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Abstract. Recently, several authors have tried to extend the usual concepts of thermodynamics 
and kinetic theory in order to deal with distributions that can be non-Boltzmannian. For dissipative 
systems described by the canonical ensemble, this leads to the notion of nonlinear Fokker-Planck 
equation (T.D. Frank, Non Linear Fokker-Planck Equations, Springer, Berlin, 2005). In this paper, 
we review general properties of nonlinear mean field Fokker-Planck equations, consider the passage 
from the generalized Kramers to the generalized Smoluchowski equation in the strong friction limit, 
and provide explicit examples for Boltzmann, Tsallis and Fermi-Dirac entropies. 
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1. THE GENERALIZED KRAMERS EQUATION 

1.1. Generalized stochastic processes 

Nonlinear Fokker-Planck (NFP) equations have been the subject of recent activity 
tfE 0, HL H. Ill] • Here, we consider a generalized Kramers equation of the form 

where h{f) and g{f) are positive functions. For /?(/) = 1 and g(/) = /, Eq. ([T]) reduces 
to the familiar Kramers equation where D is the diffusion coefficient and ^ the friction 
coefficient. Usually, 4>(r) is an external potential but we can also consider the case where 
the potential is produced by the density p (r, ? ) = / /(r, \,t)d\ according to the relation 

4>(r,f) = Ju{\r-r'\)p{r',t)dr', (2) 

where M(|r — r'|) is a binary potential of interaction. The nonlinear mean field Fokker- 
Planck equation © is associated to the Ito-Langevin stochastic process 

§ = v, ^ = -^(/)v-v$+vW)R(0, (3) 
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(4) 



where R{t) is a white noise satisfying (R(r)) = and {Ri{t)Rj{t')) = 5ij5{t - 1') where 
/ = 1 , J label the coordinates of space. 



1.2. The H-theorem 



We introduce the energy 

E = ^J fv^drd\ + ^J p^dr = K + W, (5) 

where K is the kinetic energy and W is the potential energy. For an external potential, 
we have W = J p^^dr. We define the temperature by 

T = ^. (6) 

Therefore, the Einstein relation is preserved in this generalized thermodynamical frame- 
work. We also set /3 = l/r. We introduce the generalized entropic functional 



S=-JCif)drdy, (7) 
where C(/) is a convex function (C" > 0) satisfying [@]: 



C"(/)^^|f|. (8) 

Since the temperature is fixed (canonical description), the relevant thermodynamical 
potential is the generalized free energy 

F = E-TS. (9) 

The definition of the free energy (Legendre transform) is preserved in this generalized 
thermodynamical framework. A straightforward calculation then shows that [l6|]: 

F = -J^(^Dh{f)^ + ^g{f)y^ drdy. (10) 

Therefore, F < (provided that > 0). This forms an //-theorem in the canonical 
ensemble. The free energy F plays the role of a Lyapunov functional. Note that the 
NFP equation ([T]) can be written 



df df df d 

dt dr d\ d\ 



(11) 



1.3. Stationary solutions 

The steady states of Eq. ([T]) must satisfy F = leading to a vanishing current 

J^Dh{f)^ + ^gif)y = 0. (12) 



Using Eqs. Q and ([8]), we get 



C"(/)^ + /3v = 0, (13) 
ay 



which can be integrated into 



C\f) = -P 



(14) 



where A(r) is a function of the position. Since df/dt = and J = 0, the advective 
(Vlasov) term in Eq. ([1]) must also vanish leading to the condition 



v^-V4)-^ = 0. (15) 



Using 



C"(/)|{ = -/3VA, c"(/)|^ = -/3v, (16) 

we obtain (VA — V^») ■ v = which must be true for all v. This yields VA — V4> = 0, so 
that 

A(r)=$(r) + a//3, (17) 

where oc is a constant. Therefore, the stationary solutions of Eq. ([T]) are determined by 
the relation : 

C'(/) = -/3£-a, (18) 

where e = + 4»(r) is the individual energy. Since C is convex, this equation can be 
reversed to give 

f = F{l5e + a), (19) 

where F{x) = (C')^^(— x) is a decreasing function. Thus / = /(e) is a decreasing 
function of the energy. We have /(e) = -/3 /C"{f) < 0. 



1.4. Minimum of free energy 

The critical points of free energy at fixed mass are determined by the variational 
problem 

5F + Ta5M = 0, (20) 
where a is a Lagrange multiplier. We can easily establish that 

5E = J (^ + ^^ 5fdrdy, 5S = - J C'{f)dfdrd\. (21) 



Therefore, the variational principle (l20l) gives 



C'{f) = -l5£-a, (22) 

equivalent to Eq. (fTSi) . Therefore, a stationary solution of the GK equation ([1]) is a critical 
point of free energy F[f] at fixed mass M. Furthermore, it is shown in Ref. Ililll] that 
it is linearly dynamically stable if and only if it is a minimum (at least local) of F 
at fixed mass. Note that when 4> is an external potential, we have 5^F — —Td^S = 
(l/2)r / C"{f){5f)^drd\ > so that a critical point of F is always a minimum. 

2. THE GENERALIZED SMOLUCHOWSKI EQUATION 

We restrict ourselves to the case of a constant friction so that g{f) = f and h{f) — 
fC"{f). The generalized Kramers equation ([T]) then becomes 



df df df d 

dt dr d\ d\ 



^TfC"{f)^ + fy 



(23) 



In that case, ^{f) = b, and D{f) = Df[C{f)/ f]'. Let us derive the hydrodynamic 
moments of this equation Defining the density and the local velocity by 



fdy, pu = y/vJv, (24) 
and integrating Eq. (l23l) on velocity, we get the continuity equation 

^ + V.(pu)=0. (25) 

Next, multiplying Eq. (|23l) by v, integrating on the velocity and using the continuity 
equation (|25|) . we obtain the momentum equation 

dui dui\ dPjj _(94> 



where we have defined the pressure tensor 

Pu = jfmwjdy, (27) 

where w = v — u(r,?) is the relative velocity. We now derive the generalized Smolu- 
chowski (GS) equation from the generalized Kramers (GK) equation in the strong fric- 
tion limit (see Jv], Sec. 9). For ^ +oo with fixed T, the term in parenthesis in Eq. (l23l) 
must vanish at leading order 



r/C"(/)^+/v~0. (28) 



Then, we find that the out-of-equilibrium distribution function /o(r, v,f) is determined 
by 

,,2 



C'(/o) = -/3 



(29) 



where A (r, f) is a constant of integration that is determined by the density according to 

p{r,t) = jfody = p[X{r,t)]. (30) 

Note that the distribution function fo is isotropic so that the velocity u(r,f) = 0{^^^) 
and the pressure tensor Pij = p5ij + 0{^^^) where p is given by 



p(r,0 = ^y fQV^dy = p[X{r,t) 



(31) 



Eliminating X{r,t) between Eqs. (l30l) and (|3TI) . we find that the fluid is barotropic in the 
sense that p = p{p), where the equation of state is entirely determined by the generalized 
entropy C{f). Now, considering the momentum equation (|26l) in the limit ^ +°°, we 
find that 



(32) 



Inserting this relation in the continuity equation (l25l) . we obtain the generalized Smolu- 
chowski equation 



dt 



1 . 

^(Vp + pV4)) 



(33) 



This equation can also be obtained from a Chapman-Enskog expansion in powers ofl/^ 
It monotonically decreases the free energy (s]]: 

rPpip') 



F[P] 



>/2 



■dp'dr+ljptdr, 



(34) 



which can be deduced from the free energy Q by using Eq. (|29|) to express F[f] as a 
functional F[p] = F[/o] of the density |lalv|]- A direct calculation leads to the //-theorem 



F = ~ J -{Vp + p V$)2 Jr < 0. 



(35) 



Moreover the stationary solutions of the generalized Smoluchowski equation (1331) are 
critical points of the free energy F[p] at fixed mass, satisfying 5F — a6M = where a 
is a Lagrange multiplier. This yields p'{p') / p' dp' = — 4> leading to the condition of 
hydrostatic balance 



Vp + pV4) = 0. 



(36) 



After integration, we get p = p(4>) with p'(4>) < 0. This result can also be obtained by 
integrating / = /(e) on the velocity. Finally, a steady state of the GS equation (1331) is 
linearly dynamically stable iff it is a (local) minimum of F[p] at fixed mass ^\^. 



3. EXPLICIT EXAMPLES 

3.1. Isothermal systems: Boltzmann entropy 

If we consider the Boltzmann entropy 

Ssif] = - j flnfdrdy, 
we get the ordinary Kramers equation 



df df df d 

dt dr d\ d\ 



The stationary solution is the Boltzmann distribution 

where A is determined by the conservation of mass. The equation of state is t 
mal one 

p = pT. 

In the strong friction limit, we recover the ordinary Smoluchowski equation 



dt 



1 



(rvp+pv4)) 



The free energy is the Boltzmann free energy in physical space 

F[p] = T I plnp dr+^Jp^dr. 
The stationary solution is the Boltzmann distribution in physical space 

where A' = {In/^yi^A. 



3.2. Poly tropes: Tsallis entropy 

If we consider the Tsallis (^-entropy 

Sci[f] = -^j{f-f)drdy, 
we obtain the polytropic Kramers equation 



The stationary solution is the Tsallis (or polytropic) distribution 

(^-l)/3, 



f 



1 

9-1 



(46) 



where jU is determined by the conservation of mass. The index n of the poly trope is 



+ ■ 



1 



1 



(47) 



Isothermal distribution functions are recovered in the limit q^l (i.e. n — > -\-°°). We shall 
consider ^ > so that C is convex. For ^ > 1, i.e. n> d/2, the distribution has a compact 
support (case 1) since / is defined only for £ < £,„ = qii/[{q — \)^]. For e > £,„, we set 
/ = 0. For n = d/2, f is, the Heaviside function. For q < I, the distribution is defined 
for all energies (case 2). For large velocities, it behaves like / ~ v2n-d Therefore, the 
density and the pressure are finite only for n < — 1 , i.e. J/ (J + 2) < q < 1 . This fixes the 
range of allowed parameters. The equation of state is that of a poly trope m 



Kp^, 7=1 + . 



For n> d/2 (case 1) the polytropic constant is 



K 



1 



ASd2^ 



and for n < 



n + 1 
1 (case 2), we have 

1 



_ir(j/2)r(i-d/2+R) 
r(i+n) 



-l/n 



K 



n + 1 



• 4_ir(j/2)r(-n) 



-l/n 



(48) 



(49) 



(50) 



where A=[/3|^— 1|/^]^/'-^ ^).In the strong friction limit, we get the polytropic Smolu- 
chowski equation 



dt 



The free energy is the Tsallis free energy in physical space 

K 



J (P^-P) ^^+2.1 



The stationary solution is the Tsallis distribution in physical space 

7-1 



Ky 



1 

7-1 



(51) 



(52) 



(53) 



We note that a polytropic distribution with index q in phase space yields a polytropic 
distribution with index 7=1+ 2{q — l)/[2 + d{q — 1)] in physical space. In this sense, 
Tsallis distributions are stable laws. By comparing Eq. (|46|) with Eq. (|53l) or Eqs. Q and 
(|44l) with Eq. (I52l) we note that K plays the same role in physical space as the temperature 
r = 1//3 in phase space. It is sometimes called a "polytropic temperature" [@]. 



3.3. Fermions: Fermi-Dirac entropy 



If we consider the Fermi-Dirac entropy 



— In h 

^7o ^70 



In 



drd\, 



(54) 



df df 
— +v- — 
dt dr 



1 



+ /V 



we obtain the fermionic Kramers equation 

d\ d\ 

The stationary solution is the Fermi-Dirac distribution function 
where A > is determined by the conservation of mass. The Fermi-Dirac distribution 



(55) 



(56) 



function (1561) satisfies the constraint f < r\Q which is related to the Pauli exclusion 
principle in quantum mechanics. The isothermal distribution function (l39l) is recovered 
in the non-degenerate limit f <^ rjo (valid at high temperatures). On the other hand in 
the completely degenerate limit (valid at low temperatures) the distribution is a step 
function corresponding to a polytrope of index n = d/2. The distribution in physical 
space associated with the Fermi-Dirac statistics is 



pd/2 



2 



where /„ is the Fermi integral 



Init) 



-|-oo 



-dx. 



The quantum equation of state for fermions is given in parametric form by 



pd/2 



-h_,{t) 



J/31+1 



I At) 



(57) 



(58) 



(59) 



At high temperatures we recover the classical isothermal law p = pT and at low temper- 
atures we get a poly tropic equation of state p = Kp^ with 7= and K = ■j^i^;^)'^^'^ ■ 
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